Tensor form of magnetization damping 
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A tensor form of phenomenological damping is derived for 
small magnetization motions. This form reflects basic phys- 
ical relaxation processes for a general uniformly magnetized 
particle or film. Scalar Landau-Lifshitz damping is found to 
occur only for two special cases of system symmetry. 



I. INTRODUCTION 

The dynamics of magnetization M of a single- 
domain ferromagnetic particle is usually described by the 
Landau-Lifshitz (LL) equation H: 



dM 



- 7 M x H eff 



Ml 



M x (M x H efr ). 



(1) 



Here H c g is the effective field, 7 is the gyromagnetic ra- 
tio, M s — \M\ is the saturation magnetization and a is a 
dimensionless damping parameter. The first term in (|l|) 
directly follows from microscopic equations and describes 
an averaged precession of a large number of spins. The 
second term in (Q) was introduced phenomenologically 
in Ref. [Q just from a simple geometric consideration to 
describe the magnetization damping. Later Gilbert ||] 
declared that this damping term can be rewritten as a 
"dry friction" . Here we focus on the case of small damp- 
ing (a C 1) and neglect non- uniform magnetization mo- 
tions. 

More complicated forms of phenomenological damp- 
ing were proposed by Bar'yakhtar and co-authors M, ||] 
on the base of general symmetry considerations of ex- 
change and relativistic relaxation processes. According 
to these papers, for a uniform magnetization motion "the 
crystal symmetry should influence the form of the relax- 
ation terms" and, therefore the phenomenological damp- 
ing term should contain a damping tensor with several 
damping parameters ("hierarchy of dissipative terms") 
instead of one (isotropic) damping. 

However, a general problem of magnetization damping 
can not be solved using just symmetry considerations. 
The magnetization relaxation process appears as a result 
of microscopic interactions of spins with each other and 
with phonons, conduction electrons and so on. In other 
words, a direct connection with microscopic physics must 
be found for the damping terms. The microscopic deriva- 
tion of damping has been performed, for example, in the 



case of "valence-exchange" relaxation || and magnetiza- 
tion relaxation on paramagnetic impurities A tensor 
structure of damping term follows from the results of Refs 
§and§. 

The aim of this paper is to demonstrate that the con- 
ventional relaxation term in the LL equation is inconsis- 
tent (with the exception of two special cases) with a form 
that follows from basic equations and derive a new form 
of phenomenological damping. This new tensor form of 
damping contains only one phenomenological parameter, 
which can be found experimentally. The dynamic equa- 
tion with tensor damping reflects physical relaxation pro- 
cesses for a general uniformly magnetized particle or film. 
The length of the magnetization vector is conserved and 
a rate of change of the energy proportional to the square 
of the torque. 



II. SMALL MAGNETIZATION MOTION 

Let us consider a uniformly magnetized ferromagnetic 
particle. We shall study small-amplitude motions of mag- 
netization in the vicinity of equilibrium state M =M a Z\ , 
where Z\ is the unit vector in the equilibrium direction. 
It is well-known that if we neglect the loss of energy, 
the magnetization rotation around effective field in the 
vicinity of equilibrium, in general, is elliptical. If xi and 
yi are the unit vectors corresponding to principal ellipse 
axes, the magnetic energy E can be written in the form: 



2M„ 



2NL 



(2) 



Here V is the particle volume, Hi and H2 are posi- 
tive fields, which include both microscopic and shape 
anisotropics and the external magnetic field. 

Our aim now is to derive equations of motion for M Xl 
and M Vl . We can do this by two different methods: a) 
with the help of Landau-Lifshitz Eq. ([!]) and b) with the 
help of a normal mode approach, where the relaxation is 
introduced from basic equations. 



A. Linearized LL equations 

Using (§) and <B, we can calculate the effective field 
H e ff = — 9(£/V)/9M and write down the linearized 
equations for the transverse magnetization components 

(M Z1 ~ M s ): 



1 
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di 
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(3) 



In the absence of relaxation (a = 0), the diagonal terms 
in (^) are equal to zero and the frequency of ferromag- 
netic resonance is: 



w 



(4) 



According to Eq. (||), the diagonal terms, responsi- 
ble for relaxation, in general, are different (Hi ^ H 2 ). 
The damping coefficients are equal only in two special 
cases when Hi = H 2 : 1) the case of spherical symme- 
try and 2) the case of uniaxial symmetry when the ex- 
ternal magnetic field and equilibrium magnetization are 
oriented along the easy axis. 



B. Normal mode approach 

In this approach it is convenient to introduce the clas- 
sical spin S = — 'M.V/Ti'y. Thus the energy (0) becomes 



(5) 



We shall describe small oscillations of the magnetization 
in terms of complex variables a*, a which are classical 
analogs of creation and annihilation operators of a har- 
monic oscillator and can be introduced by a Holstein- 
Primakoff transformation M for S 2 ~ —S: 



S Xl ~ (a* + a) 
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- S yi ~ (a* - a)V2S/2i. (6) 

The energy (|^) now can be rewritten in the form 

£/h = Aa*a + (B/2)(aa + a*a*), (7) 

where A = j(Hi + H 2 )/2 and B = -/(Hi - H 2 )/2. 
The dynamic precession equations are given by 

da/dt = -iAa- iBa* , da* fdt = iAa* + iBa. (8) 

The mixed terms in (§h can be eliminated by the linear 
canonical transformation 



a = uc + vc , a — uc + vc. 



(9) 



A + Uq 



B A-loq 



Thus we describe the precession in terms of the normal 
mode (c, c*) with energy of an harmonic oscillator 



£/fi = Ll)qC*C, 



(10) 



where u>q = \J A 2 — B 2 is the frequency of ferromagnetic 
resonance equivalent to Eq. (^) (A + B = 7 Hi and A — 
B = jH 2 ). The dynamic equations for c and c* are now 
independent: 



dc/dt = —iluqc, dc* / dt = iu>oC* . 



(11) 



In order to construct a damped motion for this os- 
cillator it is necessary to consider the interaction with 
a thermal bath. This implies an introduction of micro- 
scopic interactions with magnons, phonons, etc. and, in 
general, this problem must be solved using quantum sta- 
tistical methods. Here, for simplicity, we will not focus 
on a particular microscopic relaxation mechanism. Such 
a problem in a general form was solved by Lax ||. We 
shall use this result with some brief explanation. 

Let us consider the Hamiltonian of harmonic oscillator 
interacting with a thermal bath: 



H = hujotfb + ih(tfg - bg f ) + H TB - 



(12) 



Here w and b are the creation and annihilation bose op- 
erators, Htb is the thermal bath Hamiltonian, g and its 
hermitian conjugate g^ denote the thermal bath opera- 
tors, which describe a weak interaction with the oscil- 
lator. Lax analyzed the density matrix equation with 
the Hamiltonian ( |l2| ) and obtained dynamic equations 
for classical amplitudes (b) and (b') to second order in 
the interaction, where (...) is the thermal bath averag- 
ing. Denoting c = (b) and c* = (6^), we can write these 
equations as 



(d/dt + r/)c = —i(ujQ + Aw)c, 
(d/dt + jj)c* — i(ojQ + Aw)c*, 



where 



rj — iAtu = 



(13) 



(14) 



g(u) = exp(iuHrB /K)g exp(—iuli.TB /ft)- 



[..., ...] is the commutator, i] is the relaxation rate and Alo 
is the frequency shift due to interaction with the thermal 
bath (usually |Aw| <C u)q). The most important fact is 
that the equations ([l]) for classical complex amplitudes 
of damped harmonic oscillator are general. We can use 
these equations even if we do not know a microscopic re- 
laxation mechanism and find r\ and Alo from experiment 
as phenomenological parameters. 

From Eq.©, utilizing (g) and (|) with M = 
— tvyS/V, we derive linearized equations for M Xl and 



M, 



in ■ 



d 
dt 



M X1 \ _ 



/ -77 -7# 2 

I jHl -T] 



M x 



(15) 



According to (|15|), the damping of M Xl and M yi are iden- 
tical. Such intrinsic isotropy of the transverse damping 
components is seen in the Bloch-Bloembergen relaxation 
term (see, e.g., ||). We also see that the non-diagonal 
terms in (||) and (|l5|), as expected, coincide with each 
other. The diagonal terms, responsible for relaxation, are 
different. This means that the damping term in Landau- 
Lifshitz equation (|l]) is inconsistent with the basic physics 
of a damped harmonic oscillator (excluding two special 
cases mentioned above). 
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III. CONSTRUCTING THE DAMPING TERM 

We can generalize Eq.(jl|) to the form: 

^ = - 7 MxH eff -7^x[a-(MxH ea )]. (16) 
at M s 

Here a dimensionless damping tensor a is introduced. 
Note that the new tensor damping conserves the length 
of the magnetization vector (|M| = M s ) and gives a rate 
of change of the energy proportional to the square of the 
torque M x H c ff. 

The tensor a should contain all necessary information 
about symmetry of the system. Such information is in- 
cluded in the expression for the energy of the system and 
can be expressed as a tensor d 2 {£ /Vo)/<9M<9M. Thus, 
we can consider 



a— n 



d 2 {£/Vo) 
dMdM 



<9H 



cir 



dm 



(17) 



where k is a dimensionless parameter. Substituting (g) 
into (0), we obtain the damping tensor ( |l7j ) in the vicin- 
ity of equilibrium in the form: 



(18) 




Equation ( |16| ) with the tensor a (18) must be consistent 
with the Eqs. (|ll]) for small magnetization oscillations. 
Linearizing (16) and comparing with (|l|), we obtain 



k = njM s /uj 2 . 



(19) 



In the case of two (or more) stable stationary states in 
the vicinity of equilibrium we have, in general, different 
FMR frequencies and relaxation rates woi, r/i and W02, 
772, correspondingly 



IV. DISCUSSION 

In this paper a new, tensor form of damping is de- 
rived that reflects symmetry in the magnetic system en- 
ergy. The damping tensor appears from a general form 
of interaction of the normal modes of the magnetic sys- 
tem with a thermal bath. This leads to the identical 
relaxation of transverse magnetization components, as in 
the Bloch-Blocmbcrgen equations. Our analysis is ex- 
act for small oscillations about equilibrium p0|| , but the 
tensor form (as an isotropic damping) may also apply 
to large magnetization motions. The damping tensor is 
scaled by only one phenomenological damping parame- 
ter 77, which can be obtained from the experiment. It 
is demonstrated that the conventional damping term in 
the Landau-Lifshitz equation applies only for two cases 
of high symmetry. One of the most important applica- 
tions of an anisotropic damping is sure to be the case of 



thermal magnetization fluctuations M], |12|. The form 
of damping affects, so-called, fluctuation-dissipation re- 
lation and therefore changes the estimated noise level. 
More detailed study of anisotropic magnetization damp- 
ing will be published elsewhere fl3] . 

Some indications that the LLG equation does not agree 
well with experiment in a magnetic thin film is shown in 
Ref. fl4{| . In order to check the validity of the above 
new damping form it is necessary to study experimen- 
tally in detail the ferromagnetic resonance (frequency 
and linewidth) in anisotropic magnetic systems, e.g., 
films. The aim is to demonstrate that 77 (instead of a) is 
a primary relaxation rate which can depend on frequency 
u>, temperature T, external magnetic field and so on. 
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